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We report on the three Mathematica packages hexagon, CSectors , AMBRE. They are use- 
ful for the evaluation of one- and two-loop Feynman integrals with a dependence on several kine- 
matical scales. These integrals are typically needed for LHC and ILC applications, but also for 
higher order corrections at meson factories, hexagon is a new package for the tensor reduction 
of one-loop 5-point and 6-point functions with rank R = 3 and R = A, respectively; AMBRE is a 
tool for derivations of Mellin-Barnes representations; CSectors is an interface for the pack- 
age sector_decomposition and allows a convenient, direct evaluation of tensor Feynman 
integrals. 
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1. Introduction 

In recent years, we observe that higher energy, higher luminosity, higher precision, more mas- 
sive particles at LHC and ILC, but also at low energy meson factories^ lead to completely new 
demands on the efficient evaluation of Feynman diagrams, including: 

• need of some 3-point 4-point two-loop diagrams, including double boxes, 

• need of many «-point one-loop diagrams, 

with massive and massless particles participating, leading to complicated many-scale problems. 

A few of the approaches to answer the requests will be shortly introduced, concentrating on 
our own activities and on publicly available packages. They are devoted to the evaluation of L- 
loop n-point Feynman integrals of tensor rank R, with loop momenta ki, with E external legs with 
momenta p^, and with N internal lines with masses m, and propagators l/D,: 

(/7r^/2)i J D\\..D';' ...Dl' ' 



j{a} _ ^ • \ M /CI ...g /CL 

/ x^vi ^v, ^vjv ' vi-i; 



Di = (qj - mj) 



mj , (1.2) 



where we call d = 4 — 2£ the generic dimension and V; the index of the propagator. The numerator 
may contain a tensor structure, e.g.: 

= 1, )t«')tf', k'^'---k^r (1.3) 

There exists no general algorithm for the calculation of arbitrary Feynman integrals, but there 
exists a rising number of tools, and some of them are publicly available. In the next sections, I will 
comment on recent developments around three of them. 

The Feynman integrals may be evaluated in quite different ways. One may derive for them 
(systems of coupled) difference or differential equations, or one may seek a (minimal) basis of 
(scalar) master integrals, and solve only the latter ones, etc. Often the Feynman parameter rep- 
resentation is useful, which replaces the J-dimensional momentum integrations by an appropriate 
number of par ameter integrals. Feynman parameters are introduced by the representation: 

1 _ r(v, + . . . + v„) ^r' ,,4, 



D^'DJ^-.D;" r(Vi).,.r(v„) 7o '"'h (i:,D, + ...+j:„D 

with Nv = Vi + . . .Vn- The denominator ofl^"^ contains, after introduction of Feynman parameters 
Xi, the momentum dependent function with index-exponent A^v '■ 



[m 



2^-(v,+...+v„) ^ (^xyDi + ...+XNDN)-^' = {kiMijkj-2Qjkj + J)-^\ (1.5) 



'For LHC and ILC, the statement is evident. At meson factories, we have in mind e.g. luminosity determination 
with Bhabha scattering at two-loop accuracy. For details on prospects at high-luminosity <I>- and fi-factories, see e.g. the 
recent mini-review Ipl and references quoted therein. 
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Here M is an (ZjcL) -matrix, Q = Q{xi,pi,) an L-vector and J = J{xi,mf ,pc,jPe,)- The M,Q,J are 
linear in xi. The momentum integration is now simple. Shift the momenta k such that has no 
linear term ink, k = k+ {M^ ^ ) Q, and = kMk — QM ^Q + J. For a scalar Feynman integral e.g. 
one gets: 

with 

Ul{x) = (detM), (1.7) 
Fl{x) = (detM) ^2 = -{detM)J + QMQ. (1.8) 

For one-loop functions it is: Ui{x) =detM= 1 =I^x, and so Ui 'disappears'. Further, the construct 
F\{x) = —J + may be made bilinear in x,xy: Fi{x) = —J{Y,Xi) + = Y.MjXiXj. For tensor 
Feynman integrals the expressions are a little more involved, but they have the same structure: 
sums of rationals in the x, combined with non-integer powers of U (x) and F(x) [^, ^, Q]. 

2. AMBRE.m 

There is an elegant approach to the x-integrations. In the F(x) -function (1.8) one may change 
sums of monomials in x into products (as often as necessary), by Mellin-Barnes transformations, 
e.g.: 

i'oo 

After this, one may perform the x-integrations: 



L 



|i ' ' \ h I r(ai+ ••• +aN) 



Let us look at an example, the integral V614ml, see Figure ^ In a loop-by-loop approach, after 
the first momentum integration one gets here U = \ and a first f -function ( |l.8| ), which depends yet 
on one internal momentum k\ : 

fl = m^2 [X [2] +X [3] +X [4] ] -"2 - s X[2]X[4] - PR[kl+pl,m] X[1]X[2] 
- PR[kl+pl+p2, 0] X[2]X[3] - PR[kl-p2,m] X[1]X[4] 
- PR[kl, 0] X [3]X [4] , 

leading to a 7-dimensional MB-representation; after the second momentum integration, one has: 

f2 = m-^2 [X [2] +X [3] ] ^2 - s X[2]X[3] - s X[1]X[4] - 2s X[3]X[4], 



The naming convention follows ph 
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leading to another 4-dimensional integral. After several applications of Barnes' first lemma, an 
8-dimensional integral has to be treated.^ 

The package AMBRE .m is designed for a semi-automatic derivation of Mellin-Bames (MB) 
representations for Feynman diagrams; for details and examples of use see the webpage http : 
//prac . us . edu . pl/~gluza/ambre/. The package is also available from http : / /pro- 
jects . hepf orge . org/mbtools/. Version 1.0 is described in [Q], the last released version is 
1.2. We are releasing now version 2.0, which allows to construct MB -representations for two-loop 
tensor integrals. The package is yet restricted to the so-called loop-by-loop approach, which yields 
compact representations, but is known to potentially fail for non-planar topologies with several 
scales. An instructive example has been discussed in [|5]]. 

For one-scale problems, one may safely apply AMBRE . m to non-planar diagrams. For our 
example V614ml, one gets e.g. with the 8-dimensional MB -representation scetched above the 
following numerical output after running also MB.m (see also the webpage http: / /pro- 
jects .hep forge, org/mbtools / ), at s = —11: 

V614m (-i)^'^ = -0.0522082 i -0.17002 + 0.25606 £ + 4.67 6^ + ^(£3)_ (2.3) 

In simpler problems, MB -representations are a good starting point for analytical solutions, 
typically by summing multiple sums of residua. Let us take as an example the diagram F512m 
appearing in the one-loop NLO process gg —>■ qqg, shown in Figure [l], leading to (factorizing) two- 
fold infinite sums. The steps of evaluation follow closely [^, ||], and we reproduce here the result 
for the IR-divergent part in the infrared limit t' = t = 1^ + I, with m = 1 and t = {p\ + piY, in 
closed form: 

F512m(lR) = + + 7o, (2.4) 

' " sfi + ^ V- 2^' ^ ^ 

/_i = izf!^ [iog(_5) +log(v2A) + log(v4A)] , (2-6) 
ttM /I ^ 1\ n H-t,„)-^^ ^ l3nHvi/s)-^^ i-s)-^^ _ ^^.7) 

3 Stm\v2 Va) 12 rf2 12 V; lst,n 

Whereever necessary, s has to be replaced by + ie. The term Jq develops infrared endpoint singu- 
larities from the phase space integrations, due to the proportionahty of V2 = P2P3 and V4 = p^p^ to 
the gluon momentum pi,. 

Similarly, one gets for the QED pentagon F513m: 

F513m(lR) = — +7o, (2.8) 

-1,2^' V' «=o r^"j(2«+i) 



■'We made no attempt here to simplify the situation by any of the numerous tricks and reformulations etc. known to 
experts. 
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Figure 1 : The master integrals V 6 1 4 m 1 and F 5 1 2 m and F 5 1 3 m. 



(=1,2 



"=o )(2« + l) 



-3Si{n) + 2Si{\+2n) 



,(2.10) 



The answer is less singular in e, but more complicated. The inverse binomial sums with {n) and 
{2n + 1) in 7o are performed in [js]]. The expression quoted there for 7o contained instead of the 
factor {vi/s)^'^^ its expansion in e, thus developing terms depending on powers of ln(v,/5'). This 
was a disadvantage notations for the subsequent regularization of the phase space integrals, so the 
present result is more appropriate. 

3. CSectors.m 



For Euclidean kinematics, the integrand for the multi-dimensional jc-integrations is positive 
semi-definite. In numerical integrations, one has to separate the poles in J — 4, and in doing so one 
has to avoid overlapping singularities. A method for that is sector decomposition.'^ The intention 
is to separate singular regions in different variables from each other, as is nicely demonstrated by 



an example borrowed from [ 1 1 1 



1 ,-i 1 

'^^Jo ^■^xi+«e/e[x+(l-x)3;] 
^ dx dt 

+ 







dy 



dt 



(3.1) 



At several occasions, we used the package sector_decomposition [ [rO| ] (built on the 
C++ library GINAC [ll3|]) for cross checks and felt a lack of simple treatment of Feynman integrals 
with numerators. For that reason, the interface CSectors was written; it will be made publicly 
available soon. The syntax is similar to that of AMBRE. The program input for the evaluation of the 
integral V614ml is simple; we again choose m = \,s = —\\, and the topology may be read from 
the arguments of propagator functions PR: 



*There are quite a few recent papers on that, e.g. and nice reviews are given in | pl] , p^ ]. 
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<< CSectors . m 



Options [DoSectors] 

SetOptions [DoSectors, TempFileDelete -> False, SetStrategy -> C] 

nl = n2 = n3 = n4 = n5 = n6 = n7 = 1; 
m = 1 ; s = - 1 1 ; 

invariants = {pl^2 -> m^2, p2^2 -> m^2, pi p2 -> (s - 2 m^2)/2}; 

DoSectors [ { 1 } , 

{PR[kl,0,nl] PR[k2,0,n2] PR [kl+pl , m, n3 ] 

PR [kl+k2+pl, m, n5] PR [kl+k2-p2, m, n6] PR [k2-p2 , m, n7 ] } , 
{k2, kl}, invariants] [-4, 2] 

Here, the numerator is 1 (see the first argument {1} of DoSectors), and the output contains 
the functions U2 and F2: 

Using strategy C 

U = x3 x4+x3 x5+x4 x5 + x3 x6+x5 x6+x2 { x3+x4 + x6 ) +xl {x2 + x4+x5 + x6 ) 



F = xl x4^2+13 xl x4 x5+x4^2 x5+xl x5^2+x4 x5^2+13 xl x4 x6 
+ 2 xl x5 x6 + 13 x4 x5 x6+x5^2 x6+xl x6^2+x5 x6^2+x3-^2 (x4+x5 + x6) 
+x2 {x3'-2+x4^2 + 13 x4 x6+x6-^2+x3 (2 x4 + 13 x6))+x3 (x4 ^2+ {x5+x6 ) -^2 
+x4 (2x5+13 x6) ) 



Notice the presence of a ?7-f unction and the complexity of the /^-function (compared toU = \ and 
f 1 and f 2 in the loop-by-Ioop MB -approach) due to the non-sequential, direct performance of both 
momentum integrals at once. Both U and F are evidently positive semi-definite. The numerical 
result for the Feynman integral is: 

V614ml(-5)2'' = -0.052210 i - 0.17004 + 0.24634 e + 4.8773 £^ + ^(e^). (3.2) 



The numbers may be compared to ( |2.3| ). We obtained a third numerical result, also by sector 
decomposition, with the Mathematica package FIESTA [Q: 

V614ml(-5)2^ = -0.052208 ^ - 0.17002 + 0.24622 £ + 4.8746 + ^(£^). (3.3) 

The most accurate result can be obtained with an analytical representation based on harmonic 
poly logarithmic functions [|l^, 16] obtained by solving a system of differential equations [p^]: 

V614ml(-5)2^ = -0.0522082 i -0.170013 + 0.246253 £ + 4.87500 + ^(£3). (3.4) 

All displayed digits are accurate here. 
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4. hexagon.m 

The Mathematica package hexagon v. 1.0 (19 Sep 2008) was released quite recently [30 
Sep 2008] at http : / /www-zeuthen . desy . de/theory/research/CAS . html. It may 
be used for the reduction of one-loop tensor 5- and 6-point Feynman integrals up to ranks R = A and 
/? = 3, respectively, to scalar 3- and 4-point Feynman integrals. The latter have to be evaluated by 



some other package like LoopTools [18, 19] (case of only massive internal lines) or QCDloops 



[20] (general case); both these packages make also use of FF ||2T|]. The formalism underlying this 



reduction and a short decsription as well as numerical examples may be found elsewhere [E3, 23 



24] , |25| ] so that we may hold this write-up short here. We only mention that it does not use Feynman 



parameters, but is based on recurrence relations with dimensional shifts []26|]. In this approach, 
we have shown quite recently how to cancel explicitely and completely the inverse powers of the 
leading Gram determinant, hexagon was the first publicly available tensor reduction program for 
5- and 6-point Feynman integrals with arbitrary internal masses. Now, also GOLEM9 5 [ ^ became 
public, but in the released version it handles so far only massless internal particles. 

5. Summary 

We described new features of the packages AMBRE and hexagon and of the interface 

CSector s. 
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